Quantum dynamics and entanglement of a ID Fermi gas released from a trap 
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We investigate the entanglement properties of the nonequilibrium dynamics of one-dimensional 
nonintcracting Fermi gases released from a trap. The gas of N particles is initially in the ground state 
within hard-wall or harmonic traps, then it expands after dropping the trap. We compute the time 
dependence of the von Neumann and Renyi entanglement entropies and the particle fluctuations of 
spatial intervals around the original trap, in the limit of a large number N of particles. The results 
for these observables apply to one-dimensional gases of impenetrable bosons as well. 

We identify different dynamical regimes at small and large times, depending also on the initial 
condition, whether it is that of a hard-wall or harmonic trap. In particular, we analytically show 
that the expansion from hard-wall traps is characterized by the asymptotic small-time behavior 
S ~ (1/3) \n(l/t) of the von Neumann entanglement entropy, and the relation S w 7r 2 V/3 where V 
is the particle variance, which are analogous to the equilibrium behaviors whose leading logarithms 
are essentially determined by the corresponding conformal field theory with central charge c = 1. 

The time dependence of the entanglement entropy of extended regions during the expansion from 
harmonic traps shows the remarkable property that it can be expressed as a global time-dependent 
rescaling of the space dependence of the initial equilibrium entanglement entropy. 

PACS numbers: 03.65.Ud,05.30.Fk,67.85.-d 



I. INTRODUCTION 

The recent progress in the experimental activity in 
atomic physics, quantum optics and nanoscience has pro- 
vided a great opportunity to investigate the interplay 
between quantum and statistical behaviors in particle 
systems. The great ability in the manipulation of cold 
atoms PHH allows the realization of physical systems 
which are accurately described by theoretical models 
such as dilute atomic Fermi and Bose gases, Hubbard 
and Bose-Hubbard models, with different effective spa- 
tial dimensions from one to three, achieving through ex- 
perimental checks of the fundamental paradigma of the 
condensed matter physics. Experiments of cold atoms in 
optical lattices have provided a great opportunity to in- 
vestigate the unitary quantum evolution of closed many- 
body systems, exploiting their low dissipation rate which 
maintains phase coherence for a long time 0, 0|- I n 
this experimental context, the theoretical investigation 
of nonequilibrium dynamics in quantum many-body sys- 
tems, and the time evolution of the entanglement prop- 
erties characterizing the quantum correlations, is of great 
importance for a deep understanding of the fundamental 
issues of quantum dynamics, their possible applications, 
and new developments. 

We consider the nonequilibrium quantum dynamics of 
particle systems which are initially trapped within a lim- 
ited region of space by an external force, and then re- 
leased from the trap. Interesting cases are Fermi gases, or 
Bose gases with repulsive interactions, which are initially 
confined within hard-wall or harmonic traps. The free 
expansion of gases after the drop of the trap is routinely 
exploited in experiments to infer the properties of the 
initial quantum state of the particles within the trap |3| , 
observing the interference patterns of absorption images 
in the large-time ballistic regime. The time-dependence 



of the particle density and correlations have been inves- 
tigated in the literature, see, e.g., Refs. |, GHn|. In this 
paper we focus on the time evolution of the entangle- 
ment properties, and their relations with other quantum 
correlations. Our study should provide further informa- 
tion to the general issue of the entanglement properties in 
nonequilibrium dynamics after quantum quenches, which 
have recently attracted much interest; they have been in- 
vestigated in various models, and for various global and 
local quenching mechanisms, see, e.g., Refs. |13422| . 

The quantum correlations arising during the nonequi- 
librium many-body dynamics can be characterized by the 
expectation values of product of local one-particle opera- 
tors, such as the particle density, the one-particle and 
density correlations, etc., or by their integral over a 
space region A, such as the cumulants of the particle- 
number distribution within A [23l - [27| . Quantum cor- 
relations are also characterized by the fundamental phe- 
nomenon of entanglement, which gives rise to nontrivial 
connections between different parts of extended quantum 
systems [H 01 HIHH- A 

measure of entanglement is 
achieved by computing von Neumann (vN) or Renyi en- 
tanglement entropies of the reduced density matrix of 
a subsystem. One-particle correlations and bipartite en- 
tanglement entropies provide important and complemen- 
tary information of the quantum behavior of many-body 
systems, of their ground states (in particular in connec- 
tion with critical behavior) and of their nonequilibrium 
unitary evolutions under time variations of the Hamilto- 
nian, because they probe different features of the quan- 
tum dynamics. 

In this paper we consider a one-dimensional (ID) non- 
interacting Fermi gas of N particles whose initially zero- 
temperature state is the ground state within hard walls 
or in the presence of an external harmonic potential. 
This model has a wider application, because ID Bose 
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gases in the limit of strong short-ranged repulsive inter- 
actions can be mapped into a spinless fermion gas, see, 
e.g., Ref. [23. Indeed, ID Bose gases with repulsive two- 
particle short-ranged interactions, described by the Lieb- 
Liniger model [31j , become more and more nonideal with 
decreasing the particle density, acquiring fermion-likc 
properties, so that the ID gas of impenetrable bosons, or 
Tonks-Girardeau gas 32[, provides an effective descrip- 
tion of the low-density regime of ID bosonic gases [33j j . 
Due to their exact mapping, ID gases of impenetrable 
bosons and spinless fermions share the same quantum 
correlations related to the particle density, particle fluc- 
tuations of extended regions, and bipartite entanglement 
entropies of connected parts, even during nonequilibrium 
dynamics. ID systems are investigated experimentally, 
indeed the possibility of tuning the confining potential in 
experiments allows to vary the effective spatial geometry 
of the particle systems, realizing quasi- ID geometries of 
trapped quantum gases, see, e.g., Refs. [3J-|39|. 

We study the quantum correlations of the ID Fermi gas 
of N particles after the instantaneous drop of the trap, 
or during a change of the harmonic potential. We focus 
on their large- N limit, which turns out to be rapidly ap- 
proached with increasing N. In order to characterize the 
entanglement properties, we study the time dependence 
of the entanglement entropy of extended regions in prox- 
imity to the initial trap. We show that some regimes 
of the time evolution are characterized by asymptotic 
behaviors analogous to the equilibrium (ground-state) 
ones, whose leading logarithms are essentially determined 
by the corresponding conformed field theory with central 
charge c = 1. Moreover, we investigate the relations be- 
tween entanglement entropies and the distribution of the 
particle number within the same extended region. In the 
ground state of noninteracting Fermi gases the entan- 
glement entropies are asymptotically proportional to the 
particle variance for a large number of particles [H 5(| ■ 
Thus the particle variance may be considered as an ef- 
fective probe of entanglement in these class of systems, 
which should be more easily accessible to experiments. 
We show that this feature of the ground state of non- 
interacting Fermi gases is maintained in some regimes 
of the nonequilibrium dynamics after the gas is released 
from the trap. In the case of harmonic traps, the time 
dependence of the entanglement properties during the 
expansion of the Fermi gas shows the remarkable prop- 
erty that it can be expressed as a global time-dependent 
rescaling of the initial equilibrium space dependence. 

The paper is organized as follows. In Sec. [II] we report 
the many-body wave function describing the free expan- 
sion of a Fermi gas of N particles from hard- wall traps, 
and define the observables we consider. In Sec. IIIII we 
determine the large- N time evolution of several observ- 
ables, and in particular of the entanglement entropies 
and particle fluctuations of extended regions around the 
initial trap. In Sec. IIVI we consider alternative large- TV 
limits, keeping Nt or N 2 t fixed, to study the small-time 
behaviors, i.e. for t ~ 1/N and t ~ f/^V 2 , which are 



characterized by other scaling behaviors. Sec. [V] consid- 
ers the case of a Fermi gas which expands after only one 
wall drops instantaneously, thus expanding along only 
one direction. In Sec. I VII we study the nonequilibrium 
evolution of Fermi gases in a time-dependent confining 
harmonic potential, and in particular after the instanta- 
neous drop of the harmonic trap. Finally, in Sec. IVIll we 
summarize our main results and draw our conclusions. 



II. MANY-BODY WAVE FUNCTION AND 
OBSERVABLES 

A. Free expansion from hard-wall traps 

1. Both walls drop 

We consider a gas of N spinless noninteracting 
fcrmionic particles of mass m, within a hard-wall trap 
of size [— 1,1]. We set H = 1, m = 1 and I = 1, so that 
time is measured in unit of ml 2 /H. 

At t = the system is in its ground state, whose many- 
body wave function is 



^{x x , ...,x N ;t = 0) = 



f 



det[4>i(xj)] 



(1) 



where <j>k{x) are the N lowest eigenstates of the free one- 
particle Schrodinger problem with boundary conditions 
fc (-l) = 0fc(l) = O,i.e. 



cj>k(x) = sin -k(l + x) 



&k 



— k 2 

8 k ' 



1,2,... (2) 



The free expansion of the gas after the instantaneous 
drop of the walls is described by the time-dependent wave 
function 



^(xx, ...,x N ;t) 



I 



det[tpi(xj,t)] 



(3) 



where ?pi(x,t) are the one-particle wave functions with 
initial condition ^(x, 0) = 4>i(x), which can be written 
using the free propagator P as 



ipi(x,t) 



dyP(x,t;y,Q)<jH(y), 



(4) 



P(x2,t 2 ;xx,tx) 



y/l2lT{t 2 - tx) 



exp 



i(x 2 - xx) 2 



2{t 2 -tx) 



Note that they have a definite parity = (— 1) 
Eq. Q can also be written as 



k-l 



ix 2 /(2t)+imr/2 

1> n (x,t) = - - j= [f n (x,t) + (-l) 1+n f n (-x,t)], 



fn(x,t) 



dy exp 



2t 



»y 1 



7rn\ 
~2J 



(5) 



This integral can be expressed in terms of the comple- 
mentary error function or the Fresnel functions, see, e.g., 

Refs. dim. 
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2. Only one wall drops 

We also consider the case of a gas expanding after the 
instantaneous drop of only one of the walls. For simplic- 
ity, we consider an initial trap of size [0,1]. We again 
set 1 = 1. Then, after the instantaneous drop of the hard 
wall at I = 1, the gas expands along the positive real axis. 
In this case the one-particle eigenstates of the system at 
t = are 

<f> k (x) = \/2sin(7rfcE), e fe = — k 2 , k=l,2,... (6) 

The evolution of the wave function ipi(x,i) for t > 
requires the appropriate propagator Q which ensure the 
boundary condition ipk(0,t) = 0, i.e. 

tpk{x,t) = dyQ(x,t;y,0)<f> k (y), (7) 
Jo 

Q(x 2 ,t 2 ;xi,ti) = P(x 2 ,t 2 ;xi,ti) - P(-x 2 ,t 2 ;x 1 ,t 1 ). 

B. Observables 

The equal-time two point function, the particle density 
and its correlation function can be written in terms of the 
one-particle wave functions ip n (x,t), 

N 

C(x, y, t) = (c\x, t)c(y, t)) = ^ V>i(z, tf^y, t), (8) 

i=l 

N 

p(x,t) = (n(x,t)) = C(x,x,t) = \Mx,t)\ 2 , (9) 

i=l 

G n (x, y, t) = (n(x, t)n(y, t)) c = (10) 
-\C(x,y,t)\ 2 +S(x-y)C(x,y,t), 

where c(x, t) is the time-dependent fcrmionic annihila- 
tion operator and n(x,t) — c(x,tyc(x,t) is the particle- 
density operator. 

Other interesting measures of the quantum correlations 
are related to extended spatial regions, where we may 
consider the time-dependent distribution of the particle 
number and the entanglement with the rest of the sys- 
tem. In the following, we consider extended intervals 
A = [x\, x 2 ] with x\, x 2 of the size of the initial trap. For 
example, one may just consider the interval correspond- 
ing to the initial trap [— 1, 1]. 

We consider the particle-number operator of an ex- 
tended region A 

N A (t) = / dxn(x,t). (11) 

J A 

Its expectation value and connected correlation function, 
respectively 

r rn m 

n a = (n a ), (jvx) c = / n^(TN x *)>- ( 12 ) 

Ja 1= i i=l 



characterize the particle distribution within A. For this 
purpose, it is convenient to introduce the cumulants of 
the particle distribution, which can be defined through a 
generator function as [431 

Vi m) = (-idx) m Me iX * A )\x= . (13) 
In particular, the particle variance reads 

V A ee Vf } - (N 2 A ) - (N A ) 2 - / dxdyG n (x,y,t) (14) 

J A 

(the superscript m = 2 will be understood in the case of 
the particle variance). 

A measure of the entanglement of the extended region 
A with the rest of the system is provided by the Renyi 
entanglement entropies, defined as 

S { A a) = -^—lnTrp a A (15) 
1 — a 

where p A (t) is the time-dependent reduced density ma- 
trix of the subsystem A. For a — > 1, we recover the vN 
definition 

S A ee = -Tr p A In p A (16) 

(the superscript a — 1 will be understood in the case of 
the vN entanglement entropy). 

In nonintcracting Fermi gases the particle cumulants 
and the entanglement entropies of a subsystem A can be 
related to the two-point function C restricted within A, 
which we denote by C a ■ This fact holds also in nonequi- 
librium dynamics [13j . The particle number and cumu- 
lants within A can be derived using the relations (see e.g. 
Ref. (H) 

N A = Tr C A , (17) 
V] m) = {-id z ) m Q{z,C A )\ z=0 , (18) 
g(z,X) = Trln[l + (e ,:z -l)X]. (19) 

The vN and Renyi entanglement entropies can be also 
related to the two-point function C A (see Refs. [ID, |44| 
for applications to lattice systems). 

In noninteracting Fermi gases with N particles, the 
computation of the particle cumulants and entanglement 
entropies is much simplified by introducing and exploit- 
ing the information contained in the N x N overlap ma- 
trix of the one-particle wave functions (45| , as shown in 
the studies of the ground-state entanglement properties 
of homogenous Fermi gases [ill |47|. in the presence of im- 
purities and for quantum wires [48|], and in the presence 
of a space-dependent harmonic trapping potential [40| . 
As already anticipated in Ref. [Hj], the method can be 
extended to nonequilibrium quantum dynamics, by defin- 
ing the time-dependent overlap N x N matrix 

Kim(t)= / dzt(;^(z,t)tp m (z,t), n,m = l,...,N, 

J A 

(20) 



4 



where the integration is over the spatial region A, and 
involves the time evolutions tp n {x,t) of the lowest N en- 
ergy states of the one-particle Schrodinger problem at 
t = 0. The overlap matrix A and the restricted two point 
function Ca satisfy 



TrC*l = TrA 



Vfc G N, 



(21) 



which implies that the particle cumulants and the en- 
tanglement entropies can be computed from the eigen- 
values di of the N x N overlap matrix A, which are real 
and limited, a. L G (0,1). Therefore, the particle number 
and cumulants can be derived by replacing Ca with A in 
Eq. (fTS)) . In particular, 



Va = TrA(l — A), (22) 
= Tr[A - 3A 2 + 2A 3 ], (23) 
V { A i] =Tr[A-7A 2 + 12A 3 -6A 4 ], (24) 

etc.... The vN and Renyi entanglement entropies are ob- 
tained by [45[ 



N 



S A — ^ y s a (a n ), 



(25) 



where a n are the eigenvalues of A, and 
1 



s a (X) 



l-a 



In [X a + (1 — A)"] 



(26) 



and, in particular, 

Sl (A) = -AlnA-(l-A)ln(l-A) (27) 

for the vN entropy. We also mention that the determi- 
nant of the overlap matrix of a space region A provides 
the time-dependent probability to find all particles within 
A HI, indeed 



r N 

dctA= / IT dxi |\P(xi, ...x N \t)Y 
Ja n r_i 



(28) 



In the following we drop the subscript A in the quan- 
tities related to extended regions A. 



III. TIME DEPENDENCE FOR A LARGE 
NUMBER OF PARTICLES 

In this section we consider the expansion of the gas 
from the hard- wall trap [—1,1], when both walls drop 
instantaneously, corresponding to the many-body wave 
function reported in Sec. Ill A 11 We determine the time 
evolution of entanglement entropies and particle fluctua- 
tions for a large number of particles. 



A. The large-iV limit 

The large- N limit of the equal-time two-point function, 



c-z 1 ) 2 /(2t) 



C(x,y,t) = ^- fl I d Zl dz 2 e^- zl 

xe-^- z ^/( 2t Uk(zi)(j> k (z2), (29) 
can be obtained using the completeness relation [5(| 



^2<j) k {x)cj)k{y) = S(x-y) 



(30) 



k=l 



of the spectrum of the one-particle Hamiltonian at t = 0. 
We obtain the large- N limit 

C x (x,y,t) = S ^- X \ /t] e^-y^. (31) 
n(y-x) 

Numerical results at finite N, using Eq. ((29]) . show that 
C(x,y,t) approaches its large- TV limit C OQ (x,y,t) with 
0(1/ N) corrections for any t > 0. 

If we are only interested in the traces of integer powers 
of the restriction Ca of C within an extended region A, 
Eq. (f3"Tj) can be simplified dropping the phase, i.e. 



TrCt A = Tr( 



• A: 

'A > 



C A (x,y,t) 



= sin[(y - x)/t] 
n(y - x) 



(32) 
(33) 



Note that the above result does not depend on the par- 
ticular form of the confining potential. The only essential 
ingredient is that it confines the particles within a strictly 
finite region of space. For example, the above method to 
compute the large- TV limit fails in the case of a harmonic 
trap (in this case, after using the completeness relation 
we would end up with a diverging integral, calling for 
another approach to get the large- N limit), as we shall 
see later. 

Using Eq. (|3"Tj) . we can easily derive the large- N time 
evolution of the particle density 



Poo (x,t) — Coo {x, x, t) — — - 

nt 



(34) 



which means that in the formal large- TV limit the particle 
density is independent of the position. Of course, this 
regime is approached nonuniformly with respect to the 
spatial coordinate, but, as we shall is, it is rapidly reached 
around the central region of the size of the original trap, 
i.e. |x| < 1. 

Let is now consider an extended interval A — [x±, X2], 
with Xi of the size of the original trap. The average 
number of particles within A is just given by 



>(*) = / Poo{x,t) = 

J a irts 



(35) 
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FIG. 1: (Color online) The large-TV time dependence of 
the particle cumulants of the interval A = [—1,1], and, for 
comparison, their small-t asymptotic behaviors (|49[) . (|51[) and 

m. 



where 



U = t/A, 



Xl—X\. 



(36) 



Large- TV corrections are 0(1/N). The large- TV limit of 
the particle cumulants can be computed using Eq. (|18|) . 
For example, the large- TV particle variance is obtained by 



V 00 (t)=TrC A (l-C A ) 
1 r x 2 



dw\dw2 



ain{wi — W2)/t 



Wi - w 2 



(37) 
(38) 



Analogous, although more cumbersome, expressions can 
be derived for the higher cumulants, using Eq. (fPB)) with 
Ca replaced by Ca- 

In particular, the particle variance of the interval A = 
[—1,1], which coincides with the original trap, is given 
by 



2 4 
nt n^t z 



tSi(4/t) 



(39) 



t 2 (l + 1e + In(4/i) - cos(4/t) - Ci(4/t)) 



where Ci and Si are the cosine and sine integral functions. 
At small times the following asymptotic expansion holds 



Voo(t) 



1 



pn(2/f) 



cos(4/t) 2 

16tt 2 
sin(4/Q 3 

16tt 2 



7b ' 
9 2 



In 2] 

fO(t 4 ) 



(40) 



-t 2 + 0(t 4 ) 



In Fig. [T] we show curves for the particle number and the 
first few cumulants for the interval A = [—1,1] corre- 
sponding to the initial trap. 



B. The small-t behavior of bipartite entanglement 
entropies 

The small-t asymptotic behavior of the large- N limit 
of the entanglement entropies can be analytically com- 
puted. For this purpose, we note that the large- N two- 
point function Ca, cf. Eq. (|33p . can be related to an 
appropriate continuum limit of the lattice two-point func- 
tion of free fermions. Let us formally discretize the space 
within A = [xi,x^\ as 



!J 



iA 



N 1 



i,j = l,...N, 



(41) 



where A = x% — X\. Then, we consider a discretized 
version of Ca, he. the N x N matrix 



N sin[(i - j)A/r] 
A ?r(i - j) 



Nt, 



(42) 



which is identical to the two-point function of lattice 
free fermions in the thermodynamic limit without bound- 
aries [l3ll44|. 



C 



lal 



sin kp(i — j) 
w(i - j) 



(43) 



[i, j are the lattice sites), replacing the Fermi scale kp 
with A/r, and apart from a normalization N/A which 
is the analogue of the inverse lattice spacing, which we 
may formally set to one. The entanglement entropies 
of extended regions can be derived from the two-point 
function only [l3j]. Using the results of Refs. [3,[5l[ and 
the above correspondences, we derive the the asymptotic 
behavior corresponding to Eq. (|41l) . which is 



S ( a ) 



where 



In TV 



In sin ■ 



A 



0(N 



-2/cn 



(44) 





_ 1 + oT 1 






6 






= In 2 + y a , 




r°° dt 


r 6 


Va 


~ Jo T 


. 1 — a 



1 



1 



asmht/a sinht/ sinht 



1 



(45) 
(46) 

-2t 



where also the 0(N~ 2 / a ) corrections are known [461 [52|. 

Note however that the r dependence of the formula 
(|44|l is not expected to be exact at finite r = Nt in the 
nonequilibrium expansion of the gas, because it assumes 
the existence of a scaling regime at fixed r where the two- 
point function is given by Eq. (1421) , which is the large- TV 
limit at fixed t. Therefore its validity is not guaranteed 
in the large-TV limit keeping r = Nt fixed. However, it 
is expected to be valid for large values of r, and to pro- 
vide the exact asymptotic behavior in the large-T limit. 
Thus, assuming a nonsingular matching of the behavior 
for r — > oo and i — >• 0, it allows us to exactly derive 
the asymptotic expansion in the small-t limit. Note that 
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the TV and r dependence in Eq. (|44|) combines to give a 
t dependence only. The above considerations imply the 
following small-i behaviors for the vN and Renyi entan- 
glement entropies of an interval A = [a: 1,2:2] 1 



S = a pn(l/t«) + ei] - 
S< a > = c a pn(l/tj) + e c 



12 



0(t 4 ) 



(a-l)r[(a-l)/(2a)] 



(47) 
(48) 



respectively for the vN and Renyi entanglement en- 
tropies, where tg = t/A. Analogous results can be 
derived for the particle cumulants, using the large- 
TV ground-state results for homogenous gases without 
boundaries [23 - We obtain 



V = \ [Hl/ts) + 
W2 — 1 + 7b + In 2 
and for the higher cumulants 

V (2k+l) = o(t 0j 
V2k 



o(t°), 



21, ! 



k > 1, 
o(i°) fc > 2, 



(49) 
(50) 

(51) 
(52) 



where v A = -0.0185104..., v 6 = 0.00808937..., etc... Note 
that, in the case of the particle variance, Eq. (14T)|) is in 
agreement with the asymptotic small-i expansion (14T))) 
obtained by the exact large- TV expression (|39|) . In Fig. [I] 
we compare the large- TV curves with the above small-i 
asymptotic expansions. 



C. The large-time behavior of entanglement 
entropies and particle fluctuations 

The asymptotic large-t behavior of the observables 
considered can be derived by replacing the large-i be- 
havior of the one-particle wave functions, cf. Eq. (j4j, 



2 1 - (-1)' 



(53) 



in the many-body wave function, or in the overlap matrix 
([20]) . Note that the above approximations of ip n (x, t) are 
independent of x, corresponding to the fact that when 
x <C vpt ~ Nt the one-particle wave functions within 
the interval can be approximated by a constant. Thus 
the overlap matrix reads 



nm \Xl , X2 , t ) 



1 



$ mn 



[1 -(-!)"][! -(-l) m ], (54) 



1 These results can be derived by using the formulas of Ref. [46| for 
the asymptotic large- ./V expansion of the entanglement entropies 
in systems with periodic boundary conditions, by sending N — > 
00, e/L -> keeping Nirl/L = A/t = l/t s fixed. 
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FIG. 2: (Color online) The particle number of the interval 
[-1,1]. The line shows the large-TV limit n(t) = 2/{irt), cf. 
Eq. (55). 



where tg = t/{x% — x%). This implies that the large-i 
regime of the overlap matrix is characterized by only one 
nonzero eigenvalue a\ for any TV, given by 



«i 



irH s 



N 

E 

k=l 



k 2 



±[l + 0(N-i)] (55) 



The other eigenvalues get rapidly suppressed in the large- 
t limit. Numerical results show that 



02 

a i 



o(r 



«3 

a i 



o(t-% 



(56) 



where 02,03 are the next largest eigenvalues. 

The largest eigenvalue a\ determines the asymptotic 
behaviors of all observables such as the particle number, 
particle fluctuations and entanglement entropies. We ob- 
tain the large-i asymptotic behaviors of the particle cu- 
mulants, the Renyi and vN entanglement entropies, re- 
spectively 



y(m 



a 
a- 1 



1 

ntg ' 

Ol 



(a - 1) irt s ' 



S « Oi (1 — In Oi) « [In tg + 1 + In n] . 

ntg 



(57) 
(58) 
(59) 



D. Numerical results for the interval [—1,1] at 
finite TV 

We want to check the convergence of the particle fluc- 
tuations and entanglement entropies to the large- TV be- 
haviors derived in the preceding subsections. For this 
purpose, we compute them at finite particle number 
TV using the method based on the overlap matrix, see 
Sec. lII Bl We numerically compute its eigenvalues at fixed 
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FIG. 3: (Color online) The particle variance of the the inter- 
val [—1,1] for several values of N, compared with the exact 
(full line) and small-t asymptotic (dashed line) large- iV time 
dependence derived in Sec. Mil cf. Eqs. (|39p and (|49[1 respec- 
tively (which are hardly distinguishable in the figure). 



FIG. 5: (Color online) The vN entropy of the interval [—1, 1] 
as a function of the time t. The dashed lines show the small-t 
asymptotic behavior S = ci[ln(l/fa) + ei]. The full lines show 
the curves with the next known term, i.e. Eq. (|47[) . 
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FIG. 4: (Color online) The quartic cumulant for the inter- 
val [—1, 1] as a function of the time t, compared with the 
large- N time dependence derived in Sec. Mil and the small-i 
asymptotic result ((52)) . 



FIG. 6: (Color online) The a = 2 Renyi entropy of the 
interval [—1, 1] as a function of the time t. The dashed lines 
show the small-t asymptotic behavior = C2[ln(l/fa) +e2]. 
The full lines show the curves with the next known term, i.e. 
Eq. ([48]). 



t, and then obtain the particle cumulants and the entan- 
glement entropies through Eqs. (I22II25[) . We show results 
for the interval A = [—1,1] and several values of N up to 
N = 400. Figs. HH H and M show data for the par- 
ticle number, the second and quartic cumulants, and the 
vN and a = 2 Renyi entanglement entropies. They ap- 
pear to rapidly converge to the large- N analytical results 
derived above. 

The convergence to the large- N limit keeping fixed t is 
not uniform when t — > 0. The data at fixed N shown in 
Figs. hint at some nontrivial structures at very small 
t, which get hidden (pushed toward the t = axis) by 
the large- N limit keeping t fixed. 



IV. SCALING BEHAVIORS AT SMALL TIME 

In order to investigate whether other nontrivial scal- 
ing behaviors occur at small times, we consider large- iV 
limits keeping rescaled times N K t fixed, with n > 0. As 
we shall see, the time evolution of the observables re- 
lated to extended intervals show other two distinct scal- 
ing regimes, with respect to the rescaled time variables 

t = Nt, 9 = N 2 t. (60) 
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FIG. 7: (Color online) The average particle number of the 
interval [—1, 1] as a function of r = Nt. The data appear to 
converge toward the function (f62)l . 



A. Scaling with respect to r = Nt 

We first consider the interval A = [—1,1], correspond- 
ing to the initial trap. In Fig. [7] we show data for the 
average particle number n(t) versus the scaling variable 
r = Nt. The analysis of the data up to N ~ 10 3 leads to 
the time evolution 



^(t) = ^«F n (r), 



in the large- N limit, where 



Fn(r) 



1 — 7TT / 

2/(vrr) 



for 
for 



T < 4/-7T, 
T > 4/-7T 



(61) 



(62) 



Note that the dependence for r > 4/tt corresponds to the 
large- N limit at fixed t, cf. Eq. {35]), divided by TV. The 
function i 7 ^ (r) is approached quite rapidly in the large- N 
limit, as shown in Fig. [7] We have carefully checked it for 
some specific values of r with a precision of O(10 -6 ), in 
both regions r < 7r/4 and t > 7r/4 (in particular at r = 
1/3, 1/2, 1, 4/tt, 2), by extrapolating data up to N w 10 3 
assuming I/TV corrections (more precisely, fitting them 
to a + b/N + c/N 2 ). Note that F n (r) is continuous, but 
nonanalytic at r = 4/7r, where the second derivative is 
discontinuous. It is worth mentioning that the time t = 
4/(irN), corresponding to r — 4/ir, is the time taken by a 
particle with speed kp = ttN/2, which is the Fermi scale 
of the gas, to cross the interval of size L = 2, which is 
the size of the interval considered. 

In Figs. H] and O we show some particle cumulants 
and entanglement entropies for several values of N up 
to N s» 300. Their behaviors with increasing N clearly 
identify two regions, r < 4/tt and r > 4/7T, related to the 
two distinct behaviors (|62[) of the particle-number ratio. 
For r > 4/-7T the analysis of the data of the Renyi entan- 
glement entropies show a large- N behavior substantially 
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FIG. 8: (Color online) Data for V - (l/7r 2 )lniV (bottom) 
and V (4) (top) up to N = 300, versus r = Nt. 



consistent with Eq. (14"4"1) . Indeed, 

S {a ^c a [hxN + e a +F sM ( T )}, 
and, analogously, for the particle cumulants 
1 



V 



pbaJV + w 2 + F v (t)] , 



V^ 3 ) w 0, 

V (4) ^v 4 = -0.0185104.. 



(63) 

(64) 

(65) 
(66) 



Moreover, the numerical results are consistent with 

F sM (t) = F v (t) = F(t), (67) 

as shown in Fig. [TU] by the comparison of the data for 
S^ 2 ' /c a — In N — ei and w 2 V — \nN — W2 in the region 
r > 4/?r. For sufficiently large r, r > 2, the large- N limit 
of the data is well approximated by (see Fig. ^ 



F{t) ps lnsin(2/r) 



(68) 



which is the r dependence predicted by the arguments 
reported in Sec. MI B[ cf. Eq. (14"4"1) . for large r. How- 
ever, this equation does not appear to be exact, because 
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FIG. 9: (Color online) We show subtracted data of the vN 
and a — 2 entropy, respectively 5 — (1/3) In TV and — 
(1/4) In TV. In the case of the vN entropy we show the curves 
(l/3)[ln(2/r) + ei] (dashed line) and (l/3)[lnsin(2/r) + ei] 
(full line) for comparison. 



FIG. 11: (Color online) The vN entanglement entropy 
for some specific values of r versus 1/N, after subtracting 
(l/3)lnTV (bottom) and (l/2)lnTV (top). 
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FIG. 10: (Color online) We show S {2) /c a -\n N-e 2 andvrV- 
ln N— W2 in the region r > They suggest the convergence 
to a unique curve in the large- TV limit. The dashed line shows 
the curve lnsin(2/r). 



there are small 0(t 2 ) deviations, which increases when 
approaching the value r = 4/tt, see Fig. [TU] 

The large- N limit is less clear for r < 4/ir. The data 
of the particle variance and entanglement entropies, after 
subtracting the leading large-TV behavior of the initial 
ground state, does not appear to converge to any large- 
TV curve. The difference is also clearly observed for 
see the top Fig. [FJ which passes from a stable behavior for 
t > A/tt to a behavior characterized by large oscillations. 
Fig.fTTlshows data of S for selected values of r up to TV = 
200. They show that when we subtract the static leading 
log behavior, i.e. (1/3) In TV, the data for r > 4/ir look 
stable, while those for r < A/ir show a clear trend, see the 
bottom Fig.[TlJ On the other hand, the data for t < 4/ir 
appear much more stable when we subtract (1/2) In TV. 
This may suggest that a different scaling regime applies 
in this region, although we cannot exclude that we are 
just observing a crossover phenomenon with a very slow 
convergence toward the eventual large- TV behavior. This 
point should deserve further investigation. 
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FIG. 12: (Color online) The particle cumulants V and V (3) 
and the a = 2 Renyi entanglement entropy versus 8 = N 2 t. 



B. Very short-time scaling with respect to 9 = N 2 t 

Another nontrivial scaling behavior is observed at very 
small times, in the large- N limit keeping 9 = N 2 t fixed. 
As shown by Fig. [T^J the data provide a clear evidence 
of the scaling behaviors 



« w m {9), « a a (8) 



(69) 



in the large-A limit keeping 8 fixed. Note that the 
particle-number ratio R n = n(t)/N is one in this regime, 
because 1 — R n (t) appears to vanish in the large- N limit 
keeping 9 fixed, with 0(1/ N) corrections. Therefore, this 
behavior occurs at very small time scales when even the 
effective modes at the fermi scale kp = ttN/2 are still all 
practically confined within the trap. 



C. Particle cumulants and entanglement entropies 
of the interval [-1/2,1/2]. 

We now discuss the behavior of the particle cumulants 
and entanglement entropies of an interval corresponding 
to a half of the initial trap, i.e. [-1/2,1/2]= [-1/2,1/2]. 
The average particle number within this interval at t = 
is 



n(0) 



N 
~2 



N 
~2 

1 H 



1/2 

1-1/2 
1 

2N + 



dx 



sin[7r(JV + l/2)(l + a;)] 



4sin[7r(l +x)/2] 



0{N' 2 ) 



The particle cumulants are given by [271 ] 

V (2) (0) - ^(lnA^ + w 2 ) + 0(A- 1 ) 

y(2«=+D(0) = o(N°), 
V^(0)^v 2k + o(N°) 1 k>2, 



(70) 



(71) 

(72) 
(73) 
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FIG. 13: (Color online) Particle number for the interval 
[—1/2,1/2]. The dashed lines show the small-r and large-r 
behaviors of Eq. (f76)) . 



and the Renyi entanglement entropies [45| 

S^(0) = c a (\nN + e a )+0(N~ 1 / a ) (74) 

where the constants W2, v 2 k, e Q have already been de- 
fined in Eqs flU, (HO} and (p]>. 

The large- N behavior at fixed t has been already deter- 
mined in Sec. Mil Here, we focus on the large- N behavior 
at fixed r = Nt, to check how it depends on the interval 
considered. 

The ratio between the particle number n(t) within the 
interval and the total particle N rapidly approach a large- 
N scaling function of r = Nt, as shown by the data up 
to N « 200 of Fig. [H Therefore, 



r m - n{t) 

Rn(t) = — 



F„(r), 



(75) 



Moreover, the data provide a strong evidence of simple 
behaviors in the regions r < 1/7T and r > 3/tt, i.e. 



F n (r) = 



1/2 
l/(7rr) 



for 
for 



T < 1/tT, 

t > 3/tt 



(76) 



Again, the large-r behavior corresponds to the large-A^ 
time dependence n(t) = l/(irt). 

Concerning the other observables, again the data show 
the presence of different large- A" regimes. In Fig. 03] 
we show the vN entanglement entropy after subtracting 
the corresponding t — asymptotic formula (fT4"l). For 
t > 3/-7T the data provides a strong evidence of large- 
N behaviors analogous to those of Eqs. (I63H66[) . Also 
Eq. (|68|) turns out to be valid, with 



F(t) « lnsin(l/r) 



(77) 



but again this result does not appear exact, showing 
increasing, but always very small, deviations when ap- 
proaching t = 3/tt. An analogous scaling, as in Eqs. (|63I 
I66[) . is found for r < 1/7T. On the other hand, we have 
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taining 



C A (x,y,t) = Sin[ }r X l /t] - Sin ^. + ! )A] , (78) 



n(y - x) 



n(y + x) 



after dropping an irrelevant phase. In the case of an 
interval A = [0, x], the particle number reads 



n{t) =TtC a = 



1 

Tits 



Si(2/i 4 



where 



^ = t/a 



(79) 



(80) 



FIG. 14: (Color online) We report data for the vN entan- 
glement entropy after subtraction of the asymptotic large- TV 
t = expansion (l/3)(ln TV + ei). The dashed line shows the 
function (1/3) lnsin(l/r). 
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FIG. 15: (Color online) The large- TV time dependence of the 
particle cumulants, and, for comparison, their small-t asymp- 
totic behavior, for the interval [0, 1] when only one wall is 
dropped and the gas expands along the positive real axis. 



again an intermediate region, for 1/7T < r < 3/ir where 
the data appear to favor a different In TV behavior at 
large- TV, with S « (l/2)lniV instead of (1/3) In TV. 



V. TIME EVOLUTION IN THE CASE ONLY 
ONE WALL DROPS 

We now consider the case only one wall is instanta- 
neously dropped and the gas expands along the positive 
axis, which is described by the many-body wave function 
reported in Sec. Ill A 21 We proceed as in Sec. IIII Al The 
large- TV two-point function (taken at hxed t) can be again 
derived using the completeness relation, cf. Eq. (|30[) , ob- 



Results for the large- TV limit of the other observables can 
be obtained by replacing <CA{x,y,t) in the expansion of 
Eq. (|18[) . Fig. [151 shows the large- TV limit of the particle 
number, and the particle variance for the interval [0, 1], 
corresponding to the initial trap. 

In order to compute the entanglement entropies, we 
note that the large- TV two-point function (175)) appears as 
the continuum limit of the two-point function of free lat- 
tice fermions at equilibrium in the thermodynamic limit 
with open boundary conditions, see Ref. [53]. Exploit- 
ing again this correspondence, as in Sec. IIII Al we derive 
the small-t asymptotic behaviors of the entanglement en- 
tropies of the interval A = [0, x] extracting the appropri- 
ate continuum limit from known asymptotic results for 
free fermions with open boundary conditions 0, [53| , ob- 
taining 2 



S (°0 = ^ [\n(l/t s ) + e a 

2 1 - 2 /"r[(q + l)/(2a)] 
+ (a-l)T[(a-l)/(2a)} 



In 2] 



(81) 



sm(l/ts)tl /a + 0(t 2 / a ) 



Analogously, we can derive the asymptotic small-i behav- 
ior of the particle cumulants using the results of Ref. [13] , 

V w Ml /ts) + w 2 + In 2] , (82) 

V (3) fa 0, (83) 
V (4) = « 4 /2, (84) 



etc. 



The large-t behavior can be determined as in Sec. IIII Cl 
by replacing the large-i approximation of the one-particle 
wave functions in the overlap matrix (|20|) . i.e. 



<j> n (x,t) 



2{-l) n x 



(85) 



2 These results can be derived by using the formulas of Ref. |4q | for 
the asymptotic large- N expansion of the entanglement entropies 



in systems with open boundary conditions, by sending N 
l/L -» keeping Nirl/L = A/t = l/tg fixed. 



oo, 



12 



Thus, the overlap matrix of the interval A = [0,x] is 
4x 3 (-1)™ (-l) m 



^■nm (^) 



37r 3 t 3 n 



(86) 



Again, this form of the overlap matrix has only one 
nonzero eigenvalue 



4x 3 A 1 2x 3 

n—1 



[1 + 0(1/JV)] (87) 



We thus obtain the large-i behaviors 

2 



n(t) rj ai 



97T 3 t 3 ' 



ai 



97T 3 t| ' 

a a 
ai 



a- 1 



(a - 1) 97r 3 ^ 3 : 



S ~ ai(l — lnai) 



37T 3 i 3 



lni a + 



(89) 
(90) 
(91) 



It is interesting to compare these results with those of 
the gas which freely expands along both directions, cf. 
Eqs. (|35l) and (|57ti59|) . They show different power laws, 
in particular, we note the fact that, at large times, the 
trap gets emptied much faster when only one wall drops. 



VI. FERMION GASES IN TIME-DEPENDENT 
HARMONIC TRAPS 

In this section we study the nonequilibrium evolution 
of fermion gases in a time-dependent confining harmonic 
potential, 



V{x,t) = ^K(t)x 2 , 



(92) 



starting from an equilibrium ground state configuration 
with initial trap size 



-1/2 



l — K Q 



(93) 



Particular interesting cases are the instantaneous change 
to a confining potential with different trap size If, i.e. 
n{t) = Kf for t > 0, or the complete drop of the trap, 
«(*) = for t > 0. 



where the determinant involves the one-particle wave 
functions ipj(x, t) associated with the N lowest eigensolu- 
tions of the Hamiltonian at t = 0. They are the solutions 
ij}j(x,t) of the one-particle Schrodinger equation 



.dipj{x, t) 
1 dt 



\d 2 x + V{x,t) 



(94) 



with the initial condition ipj(x, 0) = <fij{x) where <fij{x) 
are the eigensolutions of the Hamiltonian at t = 0, 
characterized by a trap size Iq, with eigenvalue Ej = 
(j + 1/2) /Iq. The solution can be expressed introducing 
a time-dependent function s(t), writing piil [56j 



ipj(x,i) = s 1 / 2 4>j(x) exp ^i-ssx' 



'■dt 1 



where 



x = x/s, 



(95) 



(96) 



(j)j(x) is the j th eigenfunction of the Schrodinger equation 
of the one-particle Hamiltonian at t = 0, thus with trap 
size /oi and the real function s(t) satisfies the nonlinear 
differential equation 



s + K,(t)s = KqS 



(97) 



with initial conditions s(0) — 1 and s(0) = 0. In App. [A] 
we report some explicit solutions of the above equation, 
for an instantaneous drop of the trap, an instantaneous 
change to a trap of different size, in the case of a linear 
time dependence of the trapping potential, and a drop of 
the trap driven by a linear time dependence. 

The above results allow us to relate the time-dependent 
nonequilibrium many-body function for t > to the equi- 
librium many-body wave function at t = 0, writing [8] 



y(xx,...,x N ;t) = s N/2 y(x u ...,x N ;0) x 



exp 



(98) 



where \&(xi, a; at; 0) is the wave function of the ground 
state of the Hamiltonian at t = 0. 



B. Self-similar time evolution of density 
correlations and spatial entanglement 



A. The time-dependent many-body wave function 

The time-dependent wave function ^ of N free spinless 
fermions, with ty(x,0) given by the ground state of the 
Hamiltonian at t = 0, can be written as I5J, k35| 



*(xi, ...,x N ;t) 



1 



det^iC^,*)], 



1. Time dependence of the particle correlations 

The equal-time two-point function can be derived using 
Eq. O, 



iV 



C(x,y,t) = ^2ipi(x,t)*ipi(y,t) = ^2i>i(x,t)*i)i(y,t), 

i=l i=l 

(99) 



13 



where 



'tj)j{x,t) — s l / 2 (j>j (x) exp ^i—ssx 2 ^j . 



This implies that 

C(x, y; t) = s~ 1 C(x, y; 0)exp 



i-^ssiy 2 - x 2 ) 



(100) 



. (ioi) 



where C(x, y; 0) is the equilibrium correlation function 
for the trap size I = lo, i.e. C(x, y;0) — C(x,y)\i—i . 

Eq. (|101j) implies that the one-body entanglement 
entropy remains unchanged during the time evolution. 
At equilibrium, the one-particle Renyi entanglement en- 
tropies cr( Q ) of the ground state of trapped free fermion 
gases increase logarithmically, indeed 



r («) = 



1 



1 - a 



\nTi Pl (x,y) a =hxN 



(102) 



where pi(x,y) = C(x,y)/N is the one-particle density 
matrix. The one-particle Renyi entropies do not change 
during the time evolution, because the time-dependent 
one-particle density matrix pi(x, y, t) = C(x, y, t) /N sat- 
isfies 



Tr Pl (x,y;t) a = Tr Pl (x, y; 0) Q . 



(103) 



The particle density p(x,t) and the current j(x,t), 
which enters the conservation law 



d t p(x;t) + d x j(x;t) = 0, 
show the simple behaviors 

p{x;t) = C(x,x,t) — —p(x;0), 
s 



N 



^2[ipi(x,t)*d x ipi(x,t) 



(104) 

(105) 
(106) 



-d x ipi(x,t)*tpi(x,t)] = -xp(x;0) 
s 

where p(x; 0) = p s (x) is the static particle density for the 
initial trap size Zo- 

The ground-state (equilibrium) particle density for a 
large number of particles is given by (setting Iq = 1) [57, 

p s (x) = N 1 / 2 [R p (0 + O(l/N)], (107) 
where C = x/N 1 / 2 , and 

R P (C) = - V2-C 2 , t<(c = V2 7 (108) 

7T 

and Rp(C) = for C > Cc — v2- By integrating the par- 
ticle density, we obtain the average particle number over 
extended intervals. For example, by integrating Eq. (|107[) 



over the symmetric interval Z = [—x, x] around the cen- 
ter of the trap, we obtain 

^M=p(0 + O(l/7V), C^x/N 1 ' 2 , (109) 



C\/2 - C 2 + 2arcsin(C/v / 2) (110) 



Therefore, using Eq. (|105[) . we obtain that the large- N 
time dependence of the average particle number with S = 
[—x, x] is given by 



N z (x,t) 
N 



■■P[C/s(t)] 



(111) 



The equal-time density-density correlation behaves as 

G n (x,y;t) = s- 2 G n (x,y;0) (112) 

where G n (x, y; 0) is the static particle density correlation 
for the initial trap size. The large- N scaling of its space 
dependence differs significantly from that of the particle 
density, indeed its large- N behavior is [59l | 



Gn(x,y) ~ NRq(N 1 ^ 2 x, N 1 ^ 2 y), 



(113) 



for x =/= y. 

A discussion of the adiabatic approximation of the uni- 
tary evolution for slow changes of the harmonic potential, 
and its limitations, can be found in Ref. 1591. 



2. Particle cumulants and entanglement entropies 

The time evolution of the bipartite entanglement 
entropy S^{x\,X2',t) and the particle cumulants 
;t) of an interval [xi,X2] can be computed 
using the method based on the overlap matrix, outlined 
in Sec. Ill Bl The time-dependent overlap matrix reads 



Aij(xi,x 2 ;t) = / dzipi(z,t)*ipj(z,t) = Aij(xi,X2;0), 

where Aij(xi, X2; 0) is the static matrix of the interval 
[5i,X2] for the initial trap size. This implies that, in the 
presence of time-dependent harmonic potential, the en- 
tanglement entropies and particle cumulants of extended 
subsystems behave as 

S {a) (x 1 ,x 2 ;t) = S^(x 1 ,X2;0), (115) 
vW(x 1 ,x 2 ;t) = vW(x 1 ,x 2 ;0). (116) 

This time dependence shows the remarkable property 
that their evolution in a time-dependent harmonic po- 
tential simply corresponds to a global rescaling of the 
space dependence within the ground state of the initial 
Hamiltonian for a trap size Iq- 

We know some asymptotic large-TV behaviors of the 
ground-state entanglement properties of Fer mi g as in the 
presence of an external harmonic potential (40j |. Setting 
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Iq = 1, the asymptotic large- N expansion of the half- 
space (i.e. of the infinite interval [—00, 0] where x = is 
the center of the trap) entanglement entropies is jioj ] 



rt(a) _ C a 
^HS - y 



\nN + \n4 + e a + 0{N- 1/a ) 



(117) 



where c a and e a are the constants reported in Eqs. (|45j) 
and (|46|) . For the half-space particle cumulants we have 



^HS 



1 

2^ 



Vr 



(2i+l) 



[ln^ + ln4 + w 2 + 0(iV" 1 )] , (118) 
o(7V°) for i > 2, (119) 



HS 

^s^^ + o^ ) fori>2, 



(120) 



where u>2 and «2i are the constants appearing in the 
Eqs. (I50H52[) . Note that the half-space particle cumulants 
and entanglement entropies remain constant during the 
time evolution. Moreover, in any free expansion from 
the harmonic trap, the entanglement of any semi-infinite 
piece [— 00, x] tends asymptotically to Sgg . 

The entanglement entropies and the particle variance 
of the symmetric interval Z = [—x, x] around the center 
of trap behave as 



S [ z\x) 



V z {x)n — [In AT- 



[lniV + e a +ln2 + / s(Q )(C)], (121) 
W2 +ln2 + ./V(C)], ( 122 ) 



where £ = x/N 1 / 2 . Moreover, the large- iV extrapolation 
of numerical (practically exact) results at fixed N turns 
out to be well described by the function 

/sw (0 = MO = lnsin(vrC/V2) + Ih(4/tt) (123) 
The higher cumulants have a much simpler large- ./V be- 



havior, i.e. V% 



(2fc+l) 



and V. 



(2k) 



V2k/2. The time de- 



pendence of these quantities during the expansion of the 
gas can be obtained by rescaling the space dependence 
of these formulas according to Eqs. (|115[) and (1116[) . 

Finally, we mention that analogous results, such as 
Eqs. ([TPS]) . (HHJ), (JUS) and (|117H122|) . apply also to a ID 
gas of impenetrable bosons in time-dependent harmonic 
traps. Moreover, these results can be straightforwardly 
extended to higher-dimensional trapped systems. One 
can easily show that the evolution of the entanglement 
entropy of connected bipartitions in a harmonic potential 
corresponds to a global rescaling of the multidimensional 
space. 



C. Large-t behavior in the case of an infinite 
expansion 

In the case on an infinite expansion of the gas, due to 
the drop of the trap, s(t) diverges for t —> 00. The results 
of the previous section imply that the entanglement en- 
tropies of any finite interval vanish in the long time limit. 



On the other hand, the entanglement entropies and par- 
ticle cumulants of any semi- infinite piece [-co, a;] tends 
asymptotically to S^(-oo, 0; 0) and V (m ^> (-co, 0; 0), 
given by Eqs. (jl21j) and (I122[) respectively. 

The large-i behavior of the entanglement entropies and 
particle cumulants of the symmetric interval Z = [—x, x] 
can be analytically inferred by observing that, since the 
time evolution is characterized by a time-dependent spa- 
tial rescaling x —t x/s(t), large time implies x/s(t) —> 0. 
Thus we should evaluate the overlap matrix for a small 
interval around the center 



dxcj) n {x)cj) m {x) » 2x0„(O)0 m (O), (124) 



which has only one nonzero eigenvalue 



N 



ii = 2xJ2M0f = 2ir / 9(7V;0) 



with 



^o) = Et-t 



N y 



1/2 \/2 



-(<-l)ir(l-i/2)a 
(-1) N 



(125) 



= (126) 



Exact calculations at fixed N and t show that 

^=0(x 2 ) 
ai 



(127) 



(128) 



where a 2 is the next largest eigenvalue. Thus, the large-t 
evolution is determined by only one eigenvalue 



a 1 (t) = ^rp(N;Q) 
s(t) 



for s(t) — > 00, 



(129) 



which determines the large-i behaviors of the observables 
which can be derived from the eigenvalues of the overlap 
matrix, such as the average particle number, the particle 
cumulants and the entanglement entropies, as shown in 

Sec. Unci 

D. Free expansion after the drop of the trap 

In this section we consider the case of a quantum 
quench after an instantaneous removal of the harmonic 
trap starting from the ground state of a harmonic po- 
tential of frequency ^/kq, for which we have the analytic 
solution 



i(t) = Vl + t 2 



(130) 



where we set ko = Iq 2 = 1- The time dependence of the 
observables considered in this paper can be easily found 
for other time dependences of the trap size, such as those 
whose scaling functions s(i) are reported in App. El by 
essentially the same steps as below. 
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FIG. 16: (Color online) The average particle number of an 
interval Q = [-6, 6] with b = y/2N 1/2 . The full line shows the 
large- iV limit (fT53)l . 



FIG. 17: (Color online) The vN entanglement entropy of 
the interval Q. The dashed line shows the function given by 
Eqs. ([135]) and (po) . 



We consider the large- N limit associated with the N- 
dependent interval 



with 



Q = [-b(N),b(N)], b{N) = V2N 1 / 2 



(131) 



around the center of the trap, so that all particles are 
initially contained within Q, at least asymptotically in 
the large- N limit. Indeed, the initial average number of 
particles within Q is given by 



Nq(0) 
N 



c 

1 

N 



(132) 



whcr 



c = 2 



2 1 / 2 |Ai'(2 1 / 2 z)| 2 - 2z|Ai(2 1 / 2 z)| 2 l (133) 



which gives c = 0.0612588.... 

In the large- N limit, the time dependence of the aver- 
age particle number within Q is obtained from Eq. (11 1 1 1) , 



N Q (t) „ N Q (t) 
N ~ N Q (0) 



p[V2/s(t)] 



(134) 



which behaves as 4/ (wt) for large i. Fig.[Tn]reports results 
obtained at finite N up to TV = 60, which show that the 
convergence to the large- N limit is quite rapid. 

The time dependence of the particle cumulants and 
entanglement entropies can be derived from the corre- 
sponding static space dependence, cf. Eqs. (|12Hll23p . by 
replacing £ with y/2/s(i). In particular, the vN entropy 
is expected to behave as 

S Q (t) = i pn JV + ei + In 2 + f Q (ir/a(t))] , (135) 



Jq(x) « In sin a: + ln(2/7r). 



(136) 



Fig.[T7]shows the vN entanglement entropy of the interval 
A up to N = 60, which appear to rapidly approach the 
above large- N time dependence (the convergence appears 
slower at small times). 

The large-time behavior is characterized by another 
scaling behavior, with respect to the time variable 



t, = t/N. 



(137) 



This is already suggested by the analysis of Sec. IVI CI 
Indeed, Eq. (|129|) gives 



ai(t) 



4 N 

7T T 



4 

Tltl 



(138) 



for the largest eigenvalue of the overlap matrix of the 
interval Q. The analysis of the numerical data at finite N 
supports it, see, e.g., the vN entanglement entropy versus 
ti shown in Fig. [TH The small-^ behavior is obtained by 
matching it with the large-t behavior given by Eq. (II 35[) , 
i.e. 

5 Q «i[ln(l/t i ) + ei+ln4] (139) 
The large-t; behavios is obtained using Eq. (|138j) . 

s Q w 4- I ln *' + 1 - H*M] (140) 

Analogous results can be derived for the Renyi entropies 
and the particle cumulants. 



3 The derivation of Eqs. (11321) and H133H uses results of Refs. 1401 , 
l59l | for the anomalous power-law large- TV scaling behavior at the 
effective boundaries of the harmonic trap. 



VII. CONCLUSIONS 

We study the nonequilibrium dynamics of a ID non- 
interacting spinless Fermi gas which is initially confined 
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FIG. 18: (Color online) Large-time scaling of the vN entan- 
glement entropy of the interval A with respect to t\ = t/N. 
The dashed and dotted lines show the asymptotic behaviors 
((139J and (fl40]) respectively. 



within a limited region of space (trap) by an external 
force, and then released from the trap. As initial condi- 
tion at t — 0, we consider a Fermi gas of N particles in 
the ground state within hard walls or in the presence of 
an external harmonic potential, as in most experimental 
realizations of cold atom systems. We study the behavior 
of the quantum correlations related to extended spatial 
regions, such as the entanglement entropy and the parti- 
cle fluctuations, after the instantaneous drop of the trap, 
or during a change of the harmonic potential. 

In order to investigate the entanglement properties 
during the expansion, we consider quantum correlations 
associated with extended regions of space in proximity to 
the initial trap, such as the vN and Renyi entanglement 
entropies, and the particle cumulants which characterize 
the distribution of the particle number within the space 
region. In order to also investigate the differences related 
to the particular quenching procedure and/or the initial 
conditions, we consider Fermi gases of N particles ini- 
tially trapped by hard walls, which freely expand after 
one or both walls drop instantaneously, and by a har- 
monic potential, which gives rise to a nonhomogenous 
initial ground state due to the space-dependence of the 
confining potential. In all cases, we focus on the behav- 
ior in the limit of a large number of particles. Different 
dynamics regimes are found during the time evolution, 
which are distinguished by focusing on the large- N limit 
keeping t fixed or keeping t times appropriate powers of 
N fixed. 

In the following we summarize the main results 
achieved by this study. 

In the case of the free expansion of a ID Fermi gas 
released from hard-wall traps located within the inter- 
val [—1,1] (we set I = 1 without loosing generality), we 
find that in large- N limit the equal-time two-point func- 
tion assumes a relatively simple form, given by Eq. (1331) . 



It turns out to be an appropriate continuum limit of the 
two-point function of a lattice free-fermion model without 
boundaries at equilibrium in the thermodynamic limit. 
This allows us to infer exact small-time asymptotic be- 
haviors in the large- N limit from corresponding asymp- 
totic expansions of bipartite entanglement entropies in 
homogeneous systems, already obtained by conformal- 
field theory and other exact methods [3 l5ll - [53l |. In 
particular, for an extended interval [xi,X2] the small-i 
behavior of the a = 1 vN and Renyi entanglement en- 
tropies of an interval are given by 

= c a pn(l/t s ) + e a ] + 0(t 2/tt ), (141) 

1 + Oi , , , , 

c a = , t 5 =t/(x 2 -xi), (142) 

where the constant e a is given by Eq. ([46)1. and also the 
0(t 2 / a ) corrections are computed, cf. Eqs. (l47l) and (|48[) . 
Note that the leading logarithmic term corresponds to 
the leading logarithmic term of the Renyi entanglement 
entropy of an interval of length I 

S [a) ^c a \n£ + b a (143) 

which is the universal behavior predicted by the confor- 
mal field theory and determined by the corresponding 
central charge c = 1 [H, • Analogous results are also 
obtained for the particle cumulants, in particular the par- 
ticle variance behaves as 

V = ± [ln(l/t s ) + w 2 + 0(t 2 )] , (144) 

see Eqs. ([55)1. ([55)1 and (gSJ for more details. Moreover, 
higher cumulants behave as lA 2fe+1 ) = o(t°) (odd cumu- 
lants) and = v 2 k + o(t°) for k > 2, cf. Eq. (|55J(. The 
large-i behaviors are also computed, and are character- 
ized by negative powers laws, see Sec. lIII CI In particular, 
the large-time behavior of the vN entropy is 

Srj — [lnt 5 + l+ln7r] (145) 

Concerning the above results a few further comments 
are in order. 

(i) They are obtained in the large- N limit keeping t 
fixed, which is not uniform when t — > 0. However, the 
analysis of numerical results at finite N, obtained using 
the method based on the overlap matrix, shows that it 
is rapidly approached with increasing N, indeed O(10 2 ) 
particles, or even less, are already sufficient to show it, 
see, e.g., Figs.UHHl 

(ii) They are obtained using only the completeness re- 
lation for the one-particle discrete spectrum of the t = 
one-particle Hamiltonian. Therefore, these results do not 
depend on the particular form of the confining potential. 
The only essential ingredient is that it confines the parti- 
cles within a strictly finite region of space. For example, 
it does not apply to a harmonic potential. 

(iii) The entanglement entropies and particle variance 
show the same asymptotic relation already found in the 
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studies of the equilibrium ground-state properties [27j, 
i.e. 

S {a) « c a ir 2 V, (146) 

which has been shown to be valid for the ground state of 
a large number of noninteracting Fermi particles, in any 
dimensions and for any subsystems, in homogeneous and 
nonhomogeneous conditions [13, Eo] • 

Analogous results are obtained in the case only one 
wall of the initial trap drops, and the gas freely expands 
along one direction only. The essential point is that the 
resulting large- N two-point function corresponds to an 
appropriate continuum limit of the two-point function 
of lattice free fermions in the thermodynamic limit with 
boundaries. Therefore, one can derive asymptotic small- 
t expansions analogous to Eqs. (|143M146|) by using the 
known results for the asymptotic expansions of the en- 
tanglement entropies and particle cumulants in homoge- 
nous systems with open boundary conditions, as shown 
in Sec. 

The convergence to the large- TV limit keeping fixed t is 
not uniform when t — > 0. This limit hides other scaling 
regimes at small times with t ~ 1/N and t ~ 1/A 2 , which 
are pushed toward the t = axis when taking the large- N 
at fixed t. They emerge when studying the large- N limit 
keeping t = Nt and 9 = N 2 t fixed, as shown in Sec. IIVI 
In particular, the large-A behavior keeping r fixed, for 
sufficiently large t (t > 4/ir for an interval equal to the 
original trap), is characterized by leading log behaviors 
analogous to those at equilibrium. On the other hand, 
the large- N scaling behavior at small r (r < 4/ir for 
an interval equal to the original trap) remains unclear, 
deserving further investigation. 

Another interesting physical case of nonequilibrium 
dynamics is that of a Fermi gas expanding from a har- 
monic trap, which is closer to the conditions of experi- 
ments with cold atoms, usually realized by trapping the 
atoms with an effective harmonic potential. The differ- 
ent initial conditions with respect to hard-wall traps give 
rise to different time dependences of the observables. 

We investigate the unitary evolution of free fermion 
gases in time-dependent harmonic traps, described by 
the potential V(x,t) — ^n(t)x 2 . We study the time de- 
pendence of one-particle observables, such as the parti- 
cle density and its correlation functions, and observables 
associated with extended space regions around the cen- 
ter of the initial trap, such as the particle fluctuations 
and the entanglement entropies. The evolution in a time 
dependent harmonic trap, starting from the equilibrium 
ground state of a given initial trap with k = k(0), show 
remarkable properties: the time dependence of al above- 
mentioned observables correspond to a global rescaling of 
the system size. For example, we prove that the Renyi en- 
tanglement entropy S z (x, t) of the interval Z — [—x,x] 
(where x — is the center of the original trap) has the 
time dependence 

S { z a \x,t) = S { z a \x/s(t),0) (147) 



where s(t) is an analytical function of the time-dependent 

potential with s(0) = 1, and 8% (x, 0) is the entangle- 
ment entropy of the interval [— x, x] of the initial equi- 
librium state. In the case of a quantum quench with 
an instantaneous removal of the harmonic potential of 
frequency t/kq, we have s(t) = yl + not 2 . The en- 
tanglement entropy (x, t) of any semi-infinite space 
B = [—oo, a:] tend asymptotically to the initial half-space 
entanglement entropy, which is the half-space entangle- 
ment entropy of the ground state in a harmonic potential 
computed in Ref. [40|, i.e. 

lim S { " ] (x, t) = S { g } (0, 0) « % (In N + e a + In 2) 

t — >oo 2 

(148) 

An analogous result applies to the particle cumulants. 
The large- A behavior of the time dependence of the en- 
tanglement entropy and particle fluctuations of the fi- 
nite interval Z — [—x,x], S z ° l \x 1 t) and V^(x,t), can 
be easily determined using Eq. (|147[) and the large-A 
space dependence of the corresponding quantity at equi- 
librium [40j, see Sec. [V] In particular we consider an 
extended interval which contains (almost) all particles 
at t = (apart from 0(1/N) corrections), and deter- 
mine the large-A time dependence of its entanglement 
entropies and particle fluctuations. We find that the 
asymptotic large-A behaviors of the entanglement en- 
tropies and particle cumulants are characterized by the 
same leading logarithms at the equilibrium, and that re- 
lation (|146p holds during the time evolution, except for 
very large times t > JV where another regime sets. 

Models of ID noninteracting spinless Fermi gases have 
a wider application, because ID Bose gases in the limit of 
strong short-ranged repulsive interactions can be mapped 
into a spinless fermion gas. The basic model to describe 
the many-body features of a boson gas confined to an 
effective ID geometry is the Lieb-Liniger model with an 
effective two-particle repulsive contact interaction [3lj |. 
The limit of infinitely strong repulsive interactions cor- 
responds to a ID gas of impenetrable bosons j3^|, the 
Tonks-Girardeau gas. ID Bose gases with repulsive two- 
particle short-ranged interactions become more and more 
nonideal with decreasing the particle density, acquiring 
fermion-like properties, so that the ID gas of impene- 
trable bosons is expected to provide an effective descrip- 
tion of the low-density regime of confined ID bosonic 
gases 33]. Therefore, due to the mapping between ID 
gases of impenetrable bosons and spinless fermions, some 
correlations in free fermion gases are identical to those of 
the hard-core boson gases, such as those related to the 
particle density, particle fluctuations of extended regions, 
and bipartite entanglement entropies of connected parts. 
Therefore, the results of this paper apply to ID repul- 
sively interacting Bose gases as well. 

A further interesting issue, worth being investigated, 
concerns the universality of the behaviors found in this 
paper, whether some of them are shared with other 
many-body systems, in particular the small-time asymp- 
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totic behaviors of the entanglement entropies and particle 
fluctuations, which resembles universal behaviors found 
at equilibrium for systems with central charge c = 1. An- 
other interesting issue concerns higher-dimensional sys- 
tems, i.e. the time-dependence of the entanglement en- 
tropies during the expansion of gas released by the two 
or three-dimensional traps. 
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Appendix A: Some analytic solutions for the 
one-particle problem in a time-dependent harmonic 

trap 



We report some solutions of the Eq. (|97l) . 
(i) In the case of an instantaneous drop of the trap, so 
that n{t) = for t > 0, 



(Al) 



(ii) Instantaneous change to a confining potential with 
trap size If, so that n(t) — lj 2 for t > 0, 



s(t) 



1 + (r 2 - 1) \sm(K /2 t/r) 



(A2) 



where r = If/la- 

(iii) Linear time dependence of the trapping poten- 
tial [59|, i.e. K,(t) = kqt and r = 1 + t, 



»(t) = [ReW(T)}- 1/2 , s(t) 



ImW(T) 



[KoReWir)} 1 / 2 ' 



(A3) 



where the complex function W(r) is the solution of the 
differential equation 



iW' = k q /2 (W 2 -t) 



(A4) 



with W(l) = 1, which can be written as a combination 
of Airy functions, 

W{T) - lK » Bi(-4/ 3 r) +C Ai(-^ 3 r)' (A5) 
«y 6 Bi(-«y 3 )-^(-«y 3 ) 
^ /6 Ai(-^ /3 )- i Ai'(-4 /3 )- 

(iv) Drop of the trap driven by a linear dependence: 

n{t) = 1 - t/td for < t < t dl (A6) 

n(t) = for t > t d . 
where we set k(0) = kq = I. In this case the scaling 
function s(t) is given by 



a(t) = [ReW(l - t/t d )]- 1/2 for < t < t d , (A7) 



and 



a(t) = y/a + b(t-t d ) + c{t-t d ) 2 for t > t d , (A8) 
with 

A + b 2 



(A9) 
(A10) 



4a ' 

a = s{t d ) 2 = 1.3067374..., 

, \ 1 1 \ 0.90519789... . 

b = 2s(t d )s'{t d ) = (All) 

t d 
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